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In response to recently renewed interests in SU(A) pure-gauge dynamics with large N, both from M/string 
duality and from finite-temperature QCD phase structure, we calculate string tensions acting between the funda¬ 
mental 4, diquark 6 and other color charges in SU(4) pure-gauge theory at temperatures below the deconfining 
phase change and above the bulk phase transition. Our results suggest 4 and 6 representations have different 
string tensions, with a ratio of (Je/a 4 ^ 1.3. We also found the deconfining phase change is not strong. 


There are renewed interests in pure-gauge dy¬ 
namics with SU(IV) gauge groups with iV > 3: 
New developments in M/string theory like Mal- 
dacena’s duality conjecture ||^ map super confor¬ 
mal field theory spectra to those in large-and 
strong-coupling pure-gauge theory ||]. Another 
line of thought is questioning how well the con¬ 
ventional large-iV expansion describes quenched 
SU(3) deconfining and related phase transitions 
or changes Q . And there is a long-standing prob¬ 
lem of how hadronic string tension is related to 
finite-temperature phase structure of hadrons. 

From a M/string-provoked viewpoint, an inter¬ 
esting quantity to look at is the ratio of string 
tensions. In SU(2) and SU(3) pure-gauge theo¬ 
ries there cannot be different string tensions: any 
irreducible representation is screened by gluons 
(in the adjoint representation) down to either sin¬ 
glet or fundamental representations: in SU(3), 
6 ® 8 = 24 © 15 © 6 0 3 . This has been con¬ 
firmed by a numerical Monte-Carlo calculation 
0. However SU(fV) gauge groups with iV > 3 
can allow different string tensions for different 
representations. For example in SU(4), the 6- 
and 10-dimensional diquark representations can¬ 
not be screened by the 15 gluon down to the 
4 fundamental: 6 ® 15 = 64 © 10 © 10 © 6 
and 10 © 15 = 70 © 64 © 10 ® 6 . So this 
theory allows at least two different string ten¬ 
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sions, one for 4 and another for 6 representa¬ 
tions. For general N these string tensions are 
classified by the A^-ality k in the center Z(A^) 
group 1^: the fundamental representation has N- 
ality k = 1, the diquark representations fc = 2, 
and the adjoint A: = 0, etc. One would ex¬ 
pect (Tfc=o = 0 and 1 < crfe^o/o'fe=i < k. In¬ 
deed naive strong coupling expansion gives the 
ratio (Tklck^i = min{fc,A^ — k} from multiple 
tiling of large Wilson loops. More sophisticated 
analysis suggests k{N — k)/{N — 1) for the non- 
super-symmetric case or sin(A;7r/A^)/ sin(7r/A^) for 
weakly broken Af=2 super-symmetry ||^ . 

From the more conventional large-A^ expan¬ 
sion, an important question yet to be resolved 
is the weakness of the deconfining phase tran¬ 
sition in the SU(3) pure-gauge thermodynamics 

. This phase transition separates the confining 
phase in which the center Z(3) symmetry is pre¬ 
served from the deconfining phase in which the 
symmetry is broken. Conventional large-anal¬ 
ysis would give the following picture [|j: 1) The 
transition is located at a temperature Td ~ 0(1). 
2) it separates the confining phase with free en¬ 
ergy F ~ 0(1) and the deconfining phase with 
F ~ 0{N^). So the phase transition would be 
a strongly first-order one with a latent heat of 
~ 0{N‘^). It should not be affected by finite Nf 
flavors of light quarks if g'^N is held fixed. Thus 
large-Af would not seem to be a good guide for 
N = 3 thermodynamics where the pure-gauge 
first-order phase transition is weak and can be 
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washed away by a few light quark flavors. On the 
other hand if large-TV were to give second-order 
phase transition, as suggested by the disappear¬ 
ance of center symmetry there |^, then it would 
still serve as a useful guide: in this case the weak 
first-order phase transition of SU(3) pure-gauge 
theory is just a peculiarity arising from the cubic 
effective interaction present only in this case. 

Also interesting is the ratio of the deconfining 
phase transition temperature and the fundamen¬ 
tal string tension. In the classic Hagedorn anal¬ 
ysis the ratio should be y^ 3 / 7 r(d — 2 ) in d space- 
time dimensions. As is well known at d=A and 
N=2 and 3 the actual ratio are smaller than this 
prediction, and so are in (i=3 0. 

In the present work we study pure-gauge SU(4) 
dynamics with single-plaquette action in the fun¬ 
damental 4 representation. We use combina¬ 
tions of Metropolis or pseudo-heatbath update 
and over-relaxation update algorithms. Our lat¬ 
tices are 4 or 6 x (4, 6 , 8 , 12, or 16)^ at more 
than a dozen values of inverse-squared coupling 
in 10.10 < (3 = 8 / 5 ^ < 11.00 with at least 1000 
sweeps and up to 8000 near the deconhnement 
phase change. We impose periodic boundary con¬ 
dition in all the four axes in order to look at finite- 
temperature systems. 

We calculate the plaquette and Polyakov lines 
in 4 (fundamental), 6 (anti-symmetric diquark), 
10 (symmetric diquark) and 15 (adjoint) rep¬ 
resentations, the deconhnement fraction derived 
from them and correlations of Polyakov lines 
in all the mentioned representations. We use the 
PPR multihit method to reduce noise. 

We conhrmed the known bulk phase transition 
at around (3 ~ 10.20 where huctuation in 
both plaquette and deconhnement fraction be¬ 
comes large. This is a phase transition that sepa¬ 
rates two diherent conhning phases with different 
correlation lengths. The bulk nature of this tran¬ 
sition is known because its position in [3 does not 
move as we change the temperature lattice size 
from Lt=4 to 6 . However with Lt=4 the change 
of correlation length drives deconhnement. So we 
choose Lt=G for the rest of our calculations. 

The average Polyakov line is the order param¬ 
eter of the center Z(4) symmetry. In the con¬ 
hning phase it vanishes, while in the deconhning 



Figure 1. Deconhnement fraction from 6 x 12^ 
lattices. Varying the threshold angle 9 does not 
change the determination of the phase change 
location: 9 = 20° (cross), 25° (diamond), 30° 
(square), hot and cold starts. 15° and 35° results 
look similar and are omitted for simplicity. 


phase it acquires a hnite value in one of the four 
Z(4) (real and imaginary) axes for 4 and 10 rep¬ 
resentations, and one of the two Z(2) (real) axes 
for 6 and 15 which are self-dual representations. 
At temperatures lower than /3=10.7, all these 
Polyakov lines cluster around the origin. Above 
that temperature they start to huctuate along the 
allowed axes and gradually deviates from the ori¬ 
gin. At temperatures above /3 = 11.0 they have 
hnite values in one of the allowed axes. The de¬ 
conhnement fraction behaves in accord: it starts 
to deviate from zero at /3=10.7, crosses the 50% 
mark at about 10.75 and reaches 100% by 11.0 
for all threshold angle between 15° and 35° (see 
Figure ||). From these we conclude that there is 
a hnite-temperature deconhning phase change at 
around /3=10.75 for Lt= 6 . On the other hand we 
do not hnd any discontinuity in the average pla¬ 
quette, and especially in the difference of space¬ 
like and timelike plaquettes, in this temperature 
range. In other words we are not hnding any ev¬ 
idence for a hrst-order phase transition for L (=6 
and volumes upto 16^. This does not exclude 
hrst-order deconhning phase transition on larger 
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Figure 2. Polyakov line correlation in the fun¬ 
damental (4, cross) and anti-symmetric diquark 
( 6 , diamond) representations at /3 = 10.65 on a 
6 x 16^ lattice with 800 configurations. The curves 
are fits to tensions of (74=0.102 and tT6=0.128. 


volumes, but it would be a weak one at best j|]. 

The correlations of dual Polyakov lines for ir¬ 
reducible representations i=4, 6 , 10 and 15 , etc, 
are related to the free energy F)(r) of a pair of in¬ 
finitely heavy dual color charges in the respective 
representations separated by the distance r = \ f]: 

Ci{r) = {Li{x)L*{x + f))s (X . 

Below the deconfining temperature, the free 
energy consists of a constant-shift, Liischer- 
Coulomb, and the string-tension terms, F)(r) = 
aiV + Ci + ai/r, so we should be able to extract the 
string tension from its long-range part. Indeed 
even with the current limited statistics we obtain 
an estimate of ( 74 = 0 . 102 ( 5 ) for the fundamental 4 
representation from our 6 x 16^ lattice at /3=10.65 
(see Figure ||). The correlation in the antisym¬ 
metric diquark 6 representation is noisier but still 
yields an estimate of ( 76 = 0 . 128 ( 28 ). So there seem 
two different string tensions. The higher 10 and 
15 representations are still too noisy to extract 


any tension estimates. 

Conclusions: For the phase structure, we con¬ 
firmed with better data that in SU(4) pure-gauge 
theory the known bulk phase transition is at 
Pb 10.2 and does not move with Lt, and for 
Lt = 6 the T 7 ^ 0 deconfining phase change is 
separated from the bulk transition and is at Pd > 
10.75. A new observation is that this T 7 ^ 0 de¬ 
confining phase change does not appear stronger 
than weakly first-order SU(3) deconfining transi¬ 
tion [pd]| . For string tensions we discovered many 
things: We see string tension signals for 4 and 6 , 
but not yet for 10 or 15 . There seem two different 
string tensions, (74 ~ 0.102(5) and aQ ~ 0.128(28) 
at /3=10.65. So their ratio is 1 < ( 75/(74 < 2. And 
combining the phase structure and tension cal¬ 
culations, Td/y^aipr = 0 ) < Td/yJcri{T To) ^ 
0.53 < y^ 3 / 7 r(d — 2). We plan to accumulate 
more statistics on larger and finer lattices, prob¬ 
ably using several motherboards of the QCDSP 
parallel supercomputer being built at the RIKEN 
BNL Research Center. 
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